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Abstract 

We study classical and quantum scattering properties in the ballistic regime of particles in 
two-dimensional chaotic billiards that are models of electron- or micro- waveguides. To this 
end we construct the purely classical counterparts of the Scattering Probability (SP) matrix 
| S„, m | 2 and Husimi distributions and specialize to the case of mixed chaotic motion. Compari- 
son between classical and quantum quantities allows us to discover the purely classical dynamical 
origin of certain general, as well as particular features that appear in the quantum description 
of the system. On the other hand, at certain values of energy the tunneling of the wavefunction 
into classically forbidden regions produces striking differences between the classical and quan- 
tum quantities. We also see the manifestation of Whispering Gallery Orbits as a self-similar 
structure in the transmission part of the classical SP matrix. A potential application of this 
phenomenom in the field of microlasers is discussed briefly. 

PACS: 05.45-a, 89.75.Kd, 42.65.Wi. 



1 Introduction 

At present, the majority of studies on the Quantum-Classical Corresponcence (QCC) of chaotic sys- 
tems concern bounded motion, for which very important results have been obtained [0, ||, |j| . Most 
works also treat the situation when chaos is fully developed. Here we would like to contribute to the 
understanding of the QCC for open systems with mixed chaotic motion, in particular wc consider 
particle motion in two-dimensional (2D) billiards whose phase space is characterized by incomplete 
(Smale) horse-shoes. As is well known, 2D billiards are popular models of ballistic mesoscopic sys- 
tems (j, H| D and microwave cavities . 

Usually and naturally, QCC is explored by means of semiclassical calculations Unfortunatelly, 
this approach poses limitations since semiclassical quantities involve the cumbersome determination 
of the action associated with the trajectories and sums which often do not converge. However, as 
Baranger and co-workers have found through their semiclassical calculations of conductance, the 
dominant contributions to the average quantum conductance are classical. This motivates us to con- 
struct purely classical counterparts of the Scattering Probability (SP) matrix | S n ,m | 2 and Husimi 
distributions for 2D open billiards. In a previous paper || we showed that knowledge of the classical 
counterpart of the SP matrix enables us to predict the global structure of the quantum SP matrix. 
Here wc further our study, specializing to the case of mixed chaotic motion, by: 1) analyzing in more 
detail the type of trajectories responsible for the formation of certain structures in the classical and 
quantum SP matrices; 2) comparing Husimi distributions with classical phase space structures; and 
3) identifying certain purely quantum effects in the behaviour of the conductance. 

This paper is organized as follows. In Section II we review the definition of the SP matrix and 
the construction of its classical counterpart. In Section III the 2D waveguide model is presented, 
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as well as its resulting quantum and classical SP matrices for a large number of modes. Section 
IV concerns the dynamical analysis of the classical SP matrix, the existence of whispering gallery 
orbits, and the structure that these create in the classical SP matrix. Some aspects of the QCC are 
studied in Section V in terms of Husimi distributions and their classical counterparts, the Transient 
Poincare Maps. The appearance of purely quantum effects in our system is analyzed in Section VI 
and in Sect. VII we make some concluding remarks. 



2 Quantum and Classical SP matrices 

Quantum SP matrix 

Quantum scattering is studied via the S matrix, S, which relates incoming to outgoing waves, 

yout = g yin^ ^ 

where the V m and V out vectors specify, respectively, waves coming into and going out of the inter- 
action region. In the case of a 2D waveguide of arbitrary shape connected to two leads, say, left (L) 
and right (R) leads, the wavefunction in the leads are of the form 

m l 

V L (x,y) = J2 K<ap(i^i) +6£exp(-iA&a;)] <fa(y) 

711 — 1 

Mr 

* R (x, y)=Y^ K *Mik*x) + b R exp(-ik R x)] <f> R (y) (2) 
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where 

is the part of the wavefunction associated with the motion perpendicular to the direction of propa- 
gation; dh ((1r) is the constant width of the left (right) lead, and Ml (Mr) is the total number of 
propagating modes supported by the left (right) lead at a given Fermi energy E. For the rest of the 
paper we shall consider dh — d,R = d, and hence Ml = Mr = M. 

Note that with this notation the S matrix and the incoming and outgoing wave vectors take the 
form 



S = ( * T ' t , ) , V m ={in \ V OUt = ( II ) , (4) 

where t ,t' , r, and r' are the transmission and reflection M x M matrices. M is the highest mode 
(the largest m beyond which the longitudinal wave vector k^ R = ^Jj^r — ^ IL ~ becomes complex). 

The symbols a L,R and b L ' R stand for the vectors a^ R and b^ R , m — 1, 2. . . M. The elements 
of the transmission (t and t') and reflection (r and r') matrices are the transmission and reflec- 
tion amplitu des give n, respectively, by t mn (E ) = (yJ k n /k m )a R /q^, t' mn (E) = (\/k n /k m )b^/b R , 
r mn (E) = (y/k n /k m )b^/a!^, and r' mn (E) = {yJk n /k m )a R /b R . The squared modulo element | i n>m | 2 
(I t'n m | 2 ) gi yes the probability amplitude for a left (right)-incoming mode m to be transmitted to 
the right (left) lead into the mode n. Similarly, | r n ,m | 2 (| r ' n . m | 2 ) i s the probability for a left 
(right)-incoming mode m to be reflected to the left (right) lead into mode n. 
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The Scattering Probability (SP) matrix is denned by the 2M by 2M matrix | S n , m | 2 , represent- 
ing the transition probability for the incoming mode m to transmit or reflect into an outgoing mode n. 

Classical SP matrix 

Given the expression for the energy in the leads 

where k m and ^ are, respectively, the longitudinal and transversal components of the total wave 
vector K with magnitude K = \f2m e E /%, the angle m between k m and K is given by 



ran 1 








— sin 1 


d\]2m e E 





(6) 



Classically a particle can enter the lead at any angle in the continuous range — \ < 9 < \, 
but since the M open modes for a given energy E are discrete, we associate a range of angles 
A0 m = m — Q m -x to each mode m. That is, we coarse-grain the classical angles. Clearly the 
classical limit is M = oo. 

We now explain the procedure to construct the purely classical counterpart of | S n ,m \ 2 ■ Consider 
a classical particle entering the cavity of the waveguide, say, from the left lead and making an angle 
9i within a range corresponding to a given mode m. The particle (ray) will generally collide with 
the walls of the cavity of the waveguide a few times before exiting the cavity either to the left or to 
the right lead, making a certain angle Of to which we can associate a mode n if Of S A6*„ . To specify 
the initial conditions for the trajectory of the i-th particle, the initial position (xi, yi) and the initial 
angle Oi must be given. In order to account for all possible types of trajectories, we take a large 
number (typically 10 5 ) of initial positions for each incoming angle Oi. By recording the number of 
particles scattered into the various ranges of associated with different outgoing modes n, we obtain 
a distribution of outgoing modes for each incoming mode m. This distribution gives the classical 
counterpart of the matrix elements | t n ,m | 2 and | r n>m | 2 of the quantum SP matrix. Similarly, to 
obtain the classical counterpart of | t' n m | 2 and | r' n m | 2 we repeat the above process but for parti- 
cles entering from the right lead. This defines the procedure for constructing the classical SP matrix. 



3 The System 

We chose the geometry of the cavity of the waveguide to be that of the "rippled" billiard, which 
consists of two hard walls: a rippled wall modulated by a periodic function y(x) = d + a£(x) 
+ L) = £(x)], and a flat wall at y — 0; d is the width of the waveguide and a the modulation 
parameter. For concretenes, the periodic function in this paper is given by £(x) = 1 — cos{2itx/ L). 
To form the waveguide, we attach two semi-infinite collinear leads of width d to a cavity whose 
length is an integer multiple of the period L. Fig. 1 shows the geometry for a two-period rippled 
cavity. 

This finite length version of the rippled billiard, a model of a quantum or electromagnetic waveg- 
uide, has been used to study certain transport manifestations of chaos in the classical and 
quantum ^| regimes. Moreover, studies of the infinitely long (periodic with no leads) rippled 
billiard, originally introduced to model beam acceleration problems Jl2[ , has also provided insight in 
the understanding of general features of periodic structures (e. g., energy band structure, structure 
of eigenfunctions, etc.) and has been utilized to explore the problem of quantum-classical correspon- 
dence of classically chaotic systems Jl3|, |l4|, . 
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Before analyzing our scattering system (cavity connected to leads) it is instructive to look at the 
dynamics of the infinitely long periodic (i.e., x mod. L) billiard. A Poincare Map (PM) of the system 
supplies us with the dynamical panorama. We choose the Poincare Surface of Section (PSS) as the 
bottom boundary, y = and the PM is generated by the iteration of the Birkhoff pair of variables 
(xj,9j), labeling the longitudinal component of the position and the angle the total momentum of 
the particle makes with the x-axis right after its j collision with the bottom wall. Previous studies 
[^[ |U)[ |ll[ [H| |TJ, have shown that the motion can be either regular, mixed, or fully chaotic 
depending on the relationship between the geometrical parameters (d, a, L). A representative picture 
of a mixed phase space (obtained with the set of parameters (d, a, L) — (1, 0.305, 5.55)) is shown in 
Fig. 2. 

In this work, we analyze the waveguide system defined by a cavity formed by a single period of the 
rippled billiard [with the same values (d, a, L) used in Fig. 2] attached to two aligned semi-infinite 
leads each of width d. The period one and period four resonance island structure of Fig. 2 is formed 
by trapped orbits within the region < x < L. In particular the central resonance is formed by 
trajectories colliding nearly perpendicular with the walls in the neighbourhood of x = L/2, the 
widest part of the cavity. Since the left (right) lead is attached to the cavity at x = (x = L), such 
island structure is also present in the waveguide system. Note that the orbits forming these islands 
are classically not accesible to scattering trajectories. 

In Figures 3 and 4 we present, respectively, the quantum and classical SP matrix for the one- 
period waveguide with M = 200 open channels. The quantum SP matrix was calculated by the 
recursive Green's function method || , and the classical SP matrix by the procedure described above 
with 10 3 ensambles (each characterized by a different 9i) of 10 5 different initial conditions. Note the 
rich structure present mainly in the transmission part (| t n>m | 2 and | t' n m | 2 ) of the SP matrix. We 
note that t = t' and r = r' , as expected from the symmetry of the system. 

The similarity between quantum and classical SP matrices is remarkable in the case of a large 
number of modes (here M — 200), suggesting that the motion belongs to the deep semiclassical 
regime. Nevertheles, we have also seen a resemblance between the quantum and classical SP matri- 
ces even for a number of modes as low as 10. Hence the calculation of the classical SP matrix as a 
tool for prediction becomes relevant in quantum regimes. 



4 Dynamical analysis of the classical SPM 

First we focus on the transmission part of the classical SP matrix which is presented in Fig. 5, this 
time as a function of angles. The first clear and interesting structure to analyze is in the region 
0.345 < 9 < 0.48, for which an enlargement is made in Fig 6. The existence of Whispering Gallery 
Orbits (WGO) is well known in cavities with concave walls. These are guided orbits along the inner 



surface of such cavities, see e.g. |16 . Note that our cavity has a concave part in L/A < x < 3L/A. 
In recent works, it has been shown that WGOs produce self-similar structures 17, lq ]. Here we can 



see that a self-similar structure also shows up in the SPM, in particular in the transmission part of 
the classical SPM (Fig. 6). 

We find that the structure of Fig. 6 is formed by trajectories that hit only on the concave part 
of the upper boundary before being transmitted: the outer part of the structure (9 > 0.465) is due 
to trajectories that hit it only once; the second generation (9 ~ [0.41,0.465]) is due to trajectories 
that hit it twice, the next generation (9 ~ [0.39,0.41]) to trajectories that hit the boundary three 
times, and so on. The point of convergence of the self-similar structure (not shown in Fig. 5 due to 
numerical lack of resolution) is 9 = 0.3325, the angle at the inflection point of the profile (x = L/A). 
A particle that would hit the upper boundary at x — L/A with 6i — 0.3325 would hit the concave 
part of the billiard an infinite number of times before leaving it at an angle 9f = 0.3325. In Fig. 7 
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whispering gallery trajectories selected from different generations of the structure of Fig. 6 are shown. 

Similarly, we can discover the types of particle motion that produce the various structures in the 
transmission part of the classical SP matrix. Fig. 8 shows typical trajectories that contribute to the 
zones marked in Fig. 5. All these orbits are very stable: small variations of the initial conditions 
(xi, yi, Oi) follow closely the same trajectory, with the same number of collisions with the upper and 
lower boundaries and with very similar final angles 9f. Noting that the marix elements | S n ,m 1 2 
gives a measure of the number of trajectories connecting a pair of angles Q{ — > Of, then for every 
high intensity structure in the classical SP matrix there corresponds a boundle of stable trajectories. 
Moreover, since all the high intensity structures present in the classical SP matrix are also present 
in its quantum counterpart we can expect the stable orbits of the underlying classical motion to play 
the dominant role in determining quantum transport properties. 



5 Transient Poincare Map (TPM) 

A complementary tool to analyze a waveguide system is the Transient Poincare Map (TPM) which 
has the same meaning as the Poincare Map (PM) but now taken for the scattering system in the 
range of interaction. The TPM (with surface of section at y = 0) is presented in Fig. 9. Fig. 9a 
shows the part of the TPM generated by trajectories whose initial conditions start in the left lead 
only and Fig. 9b shows the contribution of trajectories starting from both, right and left leads. 
Comparison of Fig. 9b with Fig. 2 shows, as expected, that the stability regions present in the PM 
for the infinite billiard become forbidden phase space regions in the TPM. 

Now we can ask about the location in the TPM of the whispering gallery orbits (WGO) . For this 
purpose we construct the TPM choosing now as surface of section the top boundary because, as ex- 
plained in Sect. IV, the WGOs are guided orbits along the concave part of the boundary. The upper 
TPM is constructed by the pairs (xj , 9j ) corresponding, respectively, to the longitudinal component 
of the position and the angle made by the total momentum with the tangent of the boundary at Xj , 
right after the j th collision with the upper wall. The location of the WGO in the TPM is shown 
(highlighted) in Fig. 10, where the numbers label the number of collisions the particles make with 
the upper boundary. A structure of this type in phace space formed by WGOs has been reported 
in p8[ . We see that this pattern evolves from top to bottom; the higher the number of collisions 
with the upper wall, the closer the angle 6' is to zero. In the limit 6' = 0, a particle will collide an 
infinite number of times along the concave part of the boundary. 

A popular tool used to explore quantum-classical correspondence in billiards is the Husimi dis- 
tribution p9|| ; it is the projection of a given quantum state onto a coherent state of minimum 
uncertainty. The Husimi distribution can be viewed as a quantum phase space probability density 
that can be directly compared with the classical phase space. See jl3| for details of the calculation 
of the Husimi distributions for the infinetly periodic rippled billiard. Here, for our waveguide system 
we calculate the Husimi distributions for each one of the m opened modes at certain Fermi energy 
and compare them with the TPM generated solely by trajectories starting with the corresponding 
A0 m . Then, using a Fermi energy that supports 20 modes we present in Figures lla-lld the TPM's 
generated by particles entering from the left lead with angles corresponding to the 2 nd , 5 , 10* , and 
20 th mode, together with the corresponding Husimi distributions. Note that Husimi distributions 
coincide with the regions of phase space with the largest density of points. This agreement suggests 
that the TPM can be a used to predict, at last qualitatively, the Husimi distributions. 
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6 Purely Quantum Features 



Recently, Ketzmerick and co-workers [pi obtained important results concerning the transport prop- 
erties of "waveguides which classically produce fully chaotic or mixed dynamics. Specifically, they 
showed that the behavior of the Landauer conductance G (given by |f J2 n S m I tn,m | 2 j where t n ^ m 
are the transmission elements of the S matrix) and the Wigner delay time can clearly distinguish 
between full and mixed chaos. See also |J. In our geometry, this distinction is exemplified in Fig. 12 
where we contrast the dimensionless conductance G, in a range of energy that supports 20 modes, 
for the cases of mixed (Fig. 12a) and globally chaotic dynamics (Fig. 12b). For the mixed case 
we used the same geometrical parameters as before [(d,a,L) = (1,0.305,5.55)] and for the globally 
chaotic case we use (d, a, L) = (1, 0.305, 2.77). The difference is clear: while for mixed dynamics the 
conductance G fluctuates strongly with sharp resonances, for global chaos it is a smooth function of 
energy. These wild fluctuations in the mixed chaos case are due to the existence of resonance and 
hierarchical states, as argued in M. 

So far we have seen that the similarity between the classical SP matrix and its quantum counter- 
part enabled us to understand various quantum features in terms of purely classical dynamics (see 
also H). We have also seen a nice agreement between the TPM and its quantum counterpart, the 
Husimi distributions. The energies for which such a good agreement occurs were chosen at random. 
However, if we select now an energy value corresponding to a sharp dip or peak in the conductance 
(mixed case) we find substancial differences between the classical and quantum quantities. As an 
example, in Fig. 13 we plot the same modes as in Fig. 11 but with an energy corresponding to 
a typical dip in the conductance, see inset of Fig. 12. Note that in contrast with Fig. 11 the 
Husimis lie predominantly on classically forbidden regions of phase space. In particular, for this 
energy they have their support on the four-period resonance. The culprit of this phenomenom is 
Heisenberg's uncertainty principle that allows the wavefunction to tunnel through the KAM barriers. 

It is also instructive to look at the wavefunctions in configuration space for the non-resonant and 
resonant cases, see Fig. 14. The plots on the left (right) of this figure correspond to the non-resonant 
(resonant) case. The non-resonant (resonant) energy is the same as that used in Fig. 11 (Fig. 13) 
for modes 2, 5 and 20. These plots reveal a striking difference. In the resonant case we notice 
a "bow-tie" pattern which shadows the classical trajectory of a particle in a four-period periodic 
orbit, corresponding to the four islands in the TPM. As mentioned above, the phase space region 
of the period four islands is not accessible to classical particle incoming from the right or left leads; 
only trajectories originated in the interior of the cavity can be trapped and form this pattern. It 
is precisely this bow-tie pattern that has been exploited experimentally in |2C[| for the construction 
of high gain microlasers with directional emission with closed resonators of high refractive index. 
We see that such a pattern can be obtained also for open cavities because quantum mechanically 
the wavefunction, due to Heisenberg's uncertainty, can penetrate into the classically forbidben ar- 
eas. Thus, as we have proposed in pl| , directional emission microlasers may also be constructed 
with open cavities. We have also computed "classical wavefunctions" for these type of cavities, con- 
structed basically in the same way as the TPM with a box-counting method to mimic intensity in 
configuration space p2| . The classical wavefuntions that we obtained look very much like the truly 
quantum wave functions for the non-resonant case except for interference patterns. Of course, the 
classical wavefunctions do not reproduce the resonance structure as this is a purely quantum effect. 



7 Concluding Remarks 

The objective of this work has been to contribute to the understanding of the classical-quantum 
correspondence as regards scattering in open billiards which serve as models of mesoscopic elec- 
tron waveguides as well as microwave cavities. Here we make a comparative analysis betweeen 
the Scattering Probability (SP) matrix and the Husimi distributions with their classical counter- 
parts, namely the classical SP matrix and the Transient Poincare Maps (TPM), respectively. As a 
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paradigm of open billiards, we use a model of a ballistic 2D waveguide formed by a rippled cavity 
attached to two colinear leads. The classical particle motion in such cavity is known to undergo the 
generic (Hamiltonian) transition to chaos jic], O, purpose of this work, 



the parameters of the cavity were chosen to produce mixed chaos (a ternary incomplete horshoe |23| 

We find a very good global similarity between classical and quantum SP matrices and, as ex- 
pected, this similarity is greater as the number of open modes in the quantum system increases. 
However, even for a moderate number of modes (~ 10 — 20) the similarity between classical and 
quantum SP matrices allows to extract important information about general as well as individual 
features of the quantum system. In particular, all high intensity patterns in the classical SP matrix, 
found to be produced by bundles of stable trajectories, appear also in the quantum SP matrix. 
Hence, knowledge of the structure of the classical SP matrix allows us to predict which incoming 
modes contribute dominantly to transmission or reflection in the quantum system || . We have also 
shown that certain self-similar structure in the transmission part of the classical SP matrix is due 
to Whispering Gallery Orbits (WGO). 

On the other hand, important differences are revealed for energy values corresponding to sharp 
peaks or dips in the conductance. For these resonance energies, the wavefunctions tunnel (thanks to 
Heisenberg's uncertainty principle) into forbidden classical phase space regions (resonance islands), 
associated with stable periodic orbits within the cavity. 

Similarly, the Husimi distributions and the TPM show an excellent agreement in the case of 
non-resonant energies. The agreement suggests the usefulness of the TPM as a tool to predict the 
support of Husimi distributions for each open mode in the system. Also, we have seen for a cer- 
tain resonant energy, where substantial differences between TPM and Husimi distributions occur, 
the wavefunctions of the system form a bow-tie pattern. We have exploited this phenomenom to 
propose the construction of microlaser resonators in open cavities [pl| . 
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FIGURE CAPTIONS 



Fig. 1 Geometry of the waveguide for a two-period rippled cavity. £(x) = 1 — cos(2nx/L). 

Fig. 2 Poincare Map for (d,a,L) = (1,0.305,5.55). The Poincare Surface of Section is located 
at y = 0. 

Fig. 3 Quantum SPM, \S n ,m\ 2 , for the one-period waveguide with (d,a,L) = (1,0.305,5.55), 
and M = 200. 

Fig. 4 Classical SPM, | 

>SVi,m| 2 1 for the one-period waveguide with (d, a, L) — (1, 0.305, 5.55), and 

M = 200. 

Fig. 5 Transmission part of the classical SPM, |i„ jm | 2 as a function of angles. 

Fig. 6 Enlargement of the transmission part of the CSM for 0.345 < 6 < 0.48. 

Fig. 7 Trajectories that contribute to the first six generations of the structure of Fig. 6 (from 
only one kick on the top boundary to six kicks before leaving the concave part of that boundary). 

Fig. 8 Typical trajectories that contribute to the zones marked in Fig. 5. 

Fig. 9 (a) Transient Poincare Map (TPM) generated by particles entering to the cavity of the 
waveguide from the left lead, (b) TPM generated by particles from both leads. In both cases the 
Poincare Surface of Section is located at y = 0. 

Fig. 10 Transient Poincare Map with Poincare Surface of Section at the top boundary and the 
location (highlighted) of the WGO that contribute to the structure of Fig. 6. 

Fig. 11 Fraction of the Transient Poincare Map with Poincare Surface of Section at y = for 
particles from the left lead (black dots) and the correspondig Husimi distribution (gray-white) for 
the (a) 2 nd , (b) 5 th , (c) 10 th , and (d) 20 th . 

Fig. 12 Dimensionless conductance G for the range of energy that supports 20 modes for the 
(a) mixed case, (d,a,L) = (1,0.305,5.55); and (b) chaotic case, (d, a, L) = (1,0.305,2.77). The 
one-period waveguide is used. 

Fig. 13 Fraction of the Transient Poincare Map with Poincare Surface of Section at y = for 
particles from the left lead (black dots) and the correspondig Husimi distribution (gray-white) for 
the (a) 2 nd , (b) 5 th , (c) 10 th , and (d) 20 th ; the resonant energy of the inset of Fig. 12 is used. 

Fig. 14 Wavefuntions for the modes (a) 2 nd , (b) 5 , and 20 th for a non-resonant (left column) 
and the resonant energy of the inset of Fig. 12 (right column). 
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